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Hamiltonian dynamical systems possessing equilibria of saddle x centre x • • • x centre stability type 
display reaction-type dynamics for energies close to the energy of such equilibria; entrance and exit from 
certain regions of the phase space is only possible via narrow bottlenecks created by the influence of 
the equilibrium points. In this paper we provide a thorough pedagogical description of the phase space 
structures that are responsible for controlling transport in these problems. Of central importance is the 
existence of a Normally Hyperbolic Invariant Manifold (NHIM), whose stable and unstable manifolds have 
sufficient dimensionality to act as separatrices, partitioning energy surfaces into regions of qualitatively 
distinct behaviour. This NHIM forms the natural (dynamical) equator of a (spherical) dividing surface 
which locally divides an energy surface into two components ('reactants' and 'products'), one on either 
side of the bottleneck. This dividing surface has all the desired properties sought for in transition state 
theory where reaction rates are computed from the flux through a dividing surface. In fact, the dividing 
surface that we construct is crossed exactly once by reactive trajectories, and not crossed by nonreactive 
trajectories, and related to these properties, minimizes the flux upon variation of the dividing surface. 

We discuss three presentations of the energy surface and the phase space structures contained in it for 
2-degree-of-freedom (DoF) systems in the threedimensional space R^, and two schematic models which 
capture many of the essential features of the dynamics for n-DoF systems. In addition, we elucidate the 
structure of the NHIM. 
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1 Introduction 

Transition state theory (TST) is one of the most important theoretical and computational approaches to 
analyzing chemical reaction dynamics, both from a qualitative and quantitative point of view. Transition 
state theory was created in the 1930's, with most of the credit being given to Eyring, Polanyi, and Wigner, 
who are referred to as the "founding trinity of TST" in Miller's important review on chemical reaction rates 
( [Miller(1998)J . Nevertheless, important contributions were also made by Evans, Farkas, Szilard, Horiuti, 
Pelzer, and Marcelin, and these are described in the discussions of the historical development of the subject 
given in [Laidler fc King(1983)[ [Pollak fc Talkner(2005)l . 

The central ideas of TST are fundamental and provide a very natural and fruitful approach to "transfor- 
mation" in many areas far removed from chemistry. For example, it has been used for certain studies in atomic 
physics Jaffe et al.(2000)| , studies of the rearrangements of clusters Komat suzaki fc Berry (1999)], stu dies of 
transport phenomena in solid state and semi-conductor physics [Jacucci et al. (1984) [ [Eckhardt (1995) , stud- 
ies of diffusion dynamics in materials Voter et al.(2QQ2)] , cosmology de Oliveira et al.(2QQ2)l , and celestial 



mechanics [Jaffe et al.(20Q2)[ [Waalkens et al.(2QQ5b)| . 

The fundamental assumptions of TST we clearly stated by Wigner, who begins by stating that he considers 
chemical reactions in a setting where the equilibrium Maxwell-Boltzmann velocity and energy distributions 
are maintained (see Mahan(1974) for a detailed discussion of this point) and for which the potential energy 
surface is known ( [Garrett (2000) J . He then gives the following assumptions from which he derives TST: 

1. the motion of the nuclei occurs on the Born-Oppenheimer potential energy surface ("electronic adia- 
baticity" of the reaction) 

2. classical mechanics adequately describes the motion of the nuclei 

3. there exists a hypersurface in phase space dividing the energy surface into a region of react ants and a 
region of products having the property that all trajectories that pass from reactants to products must 
cross this dividing surface precisely once. 

It is important to note that Wigner clearly developed his ideas in phase space^ the arena for dynamics. It 
is important to keep this in mind since a great deal of later developments occur in configuration space^ in 
which certain dynamical properties are obscured. Nevertheless, the typical mental picture for the geometry 
associated with TST is built around a saddle point of a potential energy surface for a one degree-of-freedom 
(DoF) or two DoF system, i.e. it starts with a configuration space notion of sufficiently low dimension that 
standard visualisation intuition applies. 

It follows from Wigner's assumptions that fundamental geometric structure of TST is a dividing surface 
which locally divides the energy surface into two disjoint components (which we refer to as the "bottleneck 
property") and which is free of local "recrossing" (otherwise, the flux through the dividing surface would be 
overestimated). These properties are crucial for reaction rate computations. 

The problem of how to define and construct a dividing surface with these properties was solved for 

two degrees of freedom (DoF) in the 70's and 80's by Pechukas, Pollak and others [Child fc: Pollak(1980)[ 

Pechukas(19"8T)||Pechukas fc McLafferty(1973)||Pechukas fc Pollak(1978)||Pechukas fc Pollak(1977)[|Pechukas fc Pollak(1979; 



Pollak et al.(1980)| |Pollak(1981)| [Pollak fc Pechukas (1 979) [ [Pollak fc Child(1980)l who constructed the di- 



viding surface from a periodic orbit, the so called periodic dividing surface (PODS). 

The generalization to systems with more than two degrees of freedom has posed severe problems. 
Periodic orbits lack sufficient dimensionality for the construction of dividing surfaces for systems with 
three or more degrees of freedom. Progress in constructing a dividing surface for systems with three 
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or more degrees of freedom was made with the introduction of a fundament ahy new object, a nor- 
mahy hyperbohc invariant manifold, or "NHIM" ( [ Wiggins (1 990) [ |Wiggins(1994)| ), that takes the place 
of the periodic orbit. The NHIM is the main building block (and the natural higher dimensional gener- 
alization of the periodic orbit used by Pechukas, Pollak and others) of the phase space transition state 
theory mentioned above. This has been described more fully in recent work has shown that the fun- 
damental assumptions of transition state theory can only be realized in phase space for systems with 
three or more DoF. This work is described in a series of papers [Wiggins et al.(2QQl)|,[U"zer et al.(2QQ2) 



Waalkens et al.(2QQ4a)] , [Waalkens fc Wiggins(2QQ4)j, [Waalkens et al.(2QQ4b)| , |Waalkens et al.(2QQ5a) 



Waalkens et al.(2QQ5c)|, Schubert et al.(2006) , Waalkens et al.(2008) where the fundamental dynamical 



and geometrical framework for phase space TST is developed. Morever, it is shown that new types of "phase 
space structures" are responsible for governing the dynamics. 

Chief among these structures is the NHIM. This invariant manifold is normally hyperbolic, i.e., it is of 
saddle type, having stable and unstable manifolds which are of sufficient dimensionality to act as separatrices, 
i.e. impenetrable barriers to the dynamics which serve to partition energy surfaces into regions of qualitatively 
distinct behaviour. 

In addition, the NHIM forms the natural (dynamical) equator of a (spherical) dividing surface (DS) 
which locally divides the energy surface into two components, one on either side of the bottleneck. The 
dividing surface is crossed (locally) exactly once by those trajectories which pass from one component to the 
other; it is free of local recrossing by trajectories and is of minimal directional flux, and is thus the optimal 



dividing surface sought for in variational transition state theory Waalkens & Wiggins(2004) . In a nutshell. 



the dividing surface enables a natural deflnition of the phase space regions corresponding to reactants and 
products. In chemical terms, the equator of the dividing surface, the NHIM, forms the energy surface of the 
activated complex, consisting of an oscillating supermolecule poised between the reactants and the products. 

The purpose of this paper is to develop the geometrical intuition associated with these high dimensional 
phase space structures that govern reaction dynamics by constructing, step-by-step, three representations of 
the phase space geometry for n = 2 DoF, and another model that illustrates various aspects of the phase 
space structures for an arbitrary number n of DoF. Moreover, we discuss the internal structure of the NHIM. 
Before proceeding to the develpment of these models we give a technical summary of the relevant phase space 
structures, their interperation in terms of, and implications for, reaction dynamics, and the method by which 
they are realised (the Poincare-Birkhoff normal form theory). All of the material described here can be found 
in the above references. 



2 ^^Reaction-Type" Dynamics and Phase Space Transition State 
Theory 

We consider a Hamiltonian system with n degrees of freedom, phase space coordinates {Q,P) G x and 
Hamiltonian function H. We assume that (Qo^^o) is an equilibrium point of Hamilton's equations which is 
of saddle-centre-. . .-centre stability type. By adding a constant energy term which does not change the 
dynamics we can achieve H{Qq, Pq) =0, and for simplicity of exposition, we can assume that the coordinates 
have been suitably translated so that the relevant equilibrium point (Qo, ^o) is at the origin. For much of the 
discussion below, we will consider iso-energetic geometrical structures belonging to a single positive energy 

■"^The concept of an "invariant manifold" is fundamental to the description and quantification of a wide variety of dynamical 
phenomena. However, the use of the term, and its significance, may not be so clear to someone with a non-dynamical systems 
oriented background. Therefore we give a brief explanation of the term, from which will follow its implications for trajectories in 
phase space. For our purposes we can think of the term "manifold" as being synonomous with the term "surface" . "Invariance" 
of the manifold can be characterized in two (equivalent) ways. One is that trajectories through initial conditions on the manifold 
must be contained in the manifold for their entire pasts and futures. The other way of characterizing an invariant manifold is 
by saying that the vector field giving rise to the dynamics is tangent to the manifold at every point of the manifold. It follows 
immediately from both characterizations of invariance that any trajectory with initial condition not on the invariant manifold 
can never "cross" the invariant manifold. In this way, invariant manifolds of codimension 1, i.e. of dimension one less than the 
dimension of phase space, act as separatrices. They form the "barriers to transport" or "impenetrable barriers in phase space" . 

■^We will define this more precisely shortly. However, briefly, it means that the matrix associated with the linearization 
of Hamilton's equations about this equilibrium point has two real eigenvalues of equal magnitude, with one positive and one 
negative, and n — 1 purely imaginary complex conjugate pairs of eigenvalues. We will assume that the eigenvalues satisfy a 
nonresonance condition that we will describe more fully in the following. 



3 



surface := H~^{h) for some constant h > 0. In addition, we will typically restrict our attention to 

a certain neighbourhood £, local to the equilibrium point. We will defer until later a discussion of exactly 
how this region is chosen; suffice it to say for now that the region is chosen so that a new set of coordinates 
can be constructed (the normal form coordinates) in which the Hamiltonian can be expressed (the normal 
form Hamiltonian) such that it provides an integrable nonlinear approximation to the dynamics which yields 
phase space structures to within a given desired accuracy. 

Locally, the (2n — l)-dimensional energy surface Ti{h) has the structure of S^"^~^ x R in the 2n-dimensional 
phase space. The energy surface ^{h) is split locally into two components, "reactants" and "products" , by 
a (2n — 2)-dimensional "dividing surface" that is diffeomorphic to S^"^~^ and which we therefore denote by 
S^g~^(/i). The dividing surface that we construct has the following properties:- 

• The only way that trajectories can evolve from reactants to products (and vice- versa), without leaving 
the local region £, is through S^g~^(/i). In other words, initial conditions (ICs) on this dividing surface 
specify all reacting trajectories. 

• The dividing surface that we construct is free of local recrossings; any trajectory which crosses it must 
leave the neighbourhood C before it might possibly cross again. 

• The dividing surface that we construct minimizes the flux, i.e. the directional flux through the dividing 
surface will increase upon a generic deformation of the dividing surface (see [Waalkens fc Wiggins (2004)] 
for the details). 

The fundamental phase space building block that allows the construction of a dividing surface with these 
properties is a particular Normally Hyperbolic Invariant Manifold (NHIM) which, for a fixed positive energy 
/i, will be denoted Snhim(^)- "^^^ NHIM is diffeomorphic to and forms the natural dynamical equator 

of the dividing surface: The dividing surface is split by this equator into (2n — 2)-dimensional hemispheres, 
each diffeomorphic to the open (2n — 2)-ball, We will denote these hemispheres by B'^~^'^(h) and 

^ds^b^(^) and call them the "forward reactive" and "backward reactive" hemispheres, respectively. B'^~^'^{h) 
is crossed by trajectories representing "forward" reactions (from reactants to products), while B'^^{h) is 
crossed by trajectories representing "backward" reactions (from products to reactants). 

The (2n — 3)-dimensional NHIM can be viewed as the energy surface of an (unstable) invariant subsystem 
which as mentioned above, in chemistry terminology, corresponds to the "activated complex", which as an 
oscillating "supermolecule" is located between reactants and products. 

The NHIM is of saddle stability type, having (2n — 2)-dimensional stable and unstable manifolds W^{h) 
and W^(h) that are diffeomorphic to S^^~^ x R. Being of co-dimension |^ one with respect to the energy 
surface, these invariant manifolds act as separatrices, partitioning the energy surface into "reacting" and 



"non-reacting" parts as will explain in detail in Sec. 2.2 



2.1 The normal form 

As mentioned in the previous section, reaction type dynamics are induced by equilibrium points of 
saddle X centre X . . . x centre stability type. These are equilibria for which the matrix associated with the 
linearisation of Hamilton's equations have eigenvalues which consist of a pair of real eigenvalues of equal 
magnitude and opposite sign, (+A, — A), A G R, and (n — 1) pairs of complex conjugate purely imaginary 
eigenvalues, {^iuok^ — ^^fc), c^/c ^ for /c = 2, . . . , n. 

The phase space structures near equilibria of this type exist independently of a specific coordinate sys- 
tem. However, in order to carry out specific calculations we will need to be able to express these phase space 
structures in coordinates. This is where Poincare-Birkhoff normal form theory is used. This is a well-known 
theory and has been the subject of many review papers and books, see, e.g ., [Deprit(1969) , Meyer ( 1974) [ 
Dragt fc Finn(1976)| [Arnord et al.(1988)| |Meyer(1991)| [Meyer fc Han(1992)| |Murdock(2QQ3)l . For our pur- 
poses it provides an algorithm whereby the phase space structures described in the previous section can 



^Briefly, the co-dimension of a submanifold is the dimension of the space in which the submanifold exists, minus the dimension 
of the submanifold. The significance of a submanifold being "co-dimension one" is that it is one less dimension than the space 
in which it exists. Therefore it can "divide" the space and act as a separatrix, or barrier, to transport. 
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be realised for a particular system by means of the normal form transformation which involves making a 
nonlinear symplectic change of variables, 



{q,p)=T{Q,P), (1) 

into normal form coordinates^ {q^p) which, in a local neighbourhood C of the equilibrium point, "decouples" 
the dynamics into a "reaction coordinate" and "bath modes." The coordinate transformation T is obtained 
from imposing conditions on the form of H expressed the new coordinates, {q^p)^ 

HMq,P) = H{T-\q,p)) = H{Q,P). (2) 

These conditions are chosen such that i^NF and the resulting equations of motions assume a simple form 
in which the reaction coordinate and bath modes "decouple". This decoupling is one way of understanding 
how we are able to construct the phase space structures, in the normal form coordinates, that govern the 
dynamics of reaction. 

In fact, we will assume that a (generic) non-resonance condition holds between the eigenvalues, namely 
that 

k2u;2 H h knUJn / (3) 

for all integer vectors (A:2, . . . , /c^) G Z^~^. |^ When such a condition holds, the NF procedure yields an 
explicit expression for the normalised Hamiltonian i^NF as a function of n integrals of motion: 

H^F{q,p) = K^f{I, -^2, . . . , Jn) = A/ + UJ2J2 + . . . + UJnJn + h.o.t.. (4) 

Here the higher order terms (h.o.t.) are at least of order two in the integrals (/, J2, . . . , Jn). 
The integral, /, corresponds to a "reaction coordinate" (saddle- type DoF): 

I = qiPi. (5) 

We note that there is an equivalent form of the reaction coordinate:- making the linear symplectic change of 
variables qi = {pi — q)/ a/2 and pi = {pi -\- q)/ a/2, transforms the above into the following form, which may 
be more familiar to many readers, 

I = qiPi = lipl- fi) ■ (6) 

Geometrically speaking, one can move freely between these two representations by considering the plane 
{qi^Pi) and rotating it by angle 7r/4, to give (q^p). 

The integrals J^, for /c = 2, . . . , n, correspond to "bath modes" (centre- type DoFJ^ 

Jk = ^{pl^ql) , k = 2,...,n. (7) 
In the new coordinates, Hamilton's equations have a particularly simple form: 



^We note that the inclusion of zbA in a non-resonance condition would be vacuous; one cannot have a resonance of this kind 
between a real eigenvalue, zbA, and purely imaginary eigenvalues, ±20;^, k = 2^ . . . ,n. 

^Throughout our work we use, somewhat interchangeably, terminology from both chemical reaction dynamics and dynamical 
systems theory. This is most noticable in our reference to the integrals of motion. / is the integral related to reaction, and in 
the context of dynamical systems theory it is related to hyperbolic behaviour. The term "reactive mode" might also be used 
to describe the dynamics associated with this integral. The integrals J2, • • • , Jn describe the dynamics associated with "bath 
modes" . In the context of dynamical systems theory, the dynamics associated with these integrals is referred to as "center type 
dynamics" or "center modes" . A key point here is that integrals of the motion provide us with the natural way of defining and 
describing the physical notion of a "mode". The nature of the mode is defined in the context of the specific application (i.e. 
chemical reactions) or, in the context of dynamical systems theory, through its stability properties (i.e. hyperbolic or centre). 
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Qi 
Pi 
Qk 
Pk 

for /c = 2, . . . , n, where we denote 



NF 



dpi 



NF 



dpk 
dH^F 



Ai(/,J2,. 


' ' , Jn)Ql, 


A(/,J2,.. 


',Jn)Pl, 


^k{I, J2, ' 


' ' , Jn)Pk 


^k{I, J2, ' 


• • •) Jn^Qk- 



(8) 



A(/,J2,...,J„) := ^^nf(/,J2,...,J„) 
n,{I,J„...,J^) := 9KMI,J2,...,Jn) ^ k = 2,...,n. (10) 

OJk 

The integrals provide a natural definition of the term "mode" that is appropriate in the context of reaction, 
and they are a consequence of the (local) integrability in a neighborhood of the equilibrium point of saddle- 
centre-. . .-centre stability type. Moreover, the expression of the normal form Hamiltonian in terms of the 
integrals provides us a way to partition the 'energy' between the different modes. We will provide examples 
of how this can be done in the following. 

The normal form transformation T in ([T]) can be computed in an algorithmic fashion. One can give 
explicit expression for the phase space structures discussed in the previous section in terms of the normal 
form coordinates, (q^p). This way the phase space structures can be constructed in terms on the normal form 
coordinates, and for physical interpretation, transformed back to the orriginal "physical" coordinates, 

{Q^P), by the inverse of the transformation Tj^ 

In summary, the "output" of the normal form algorithm is the following: 

• A symplectic transformation T((5, P) = {q^p)^ and its inverse T~^{q,p) = (Q, P), that relate the normal 
form coordinates (q^p) to the original "physical" coordinates {Q^P). 

• An expression for the normalized Hamiltonian: in the form, i^NF, in terms of the normal form coordi- 
nates and in the form i^NF, in terms of the integrals (/, J2, . . . , Jn). 

• Explicit expressions for the integrals of motion / and J/c, k = 2, . . . ,n, in terms of the normal form 
coordinates. 



2.2 Explicit definition of the phase space structures in the normal form coordi- 
nates 

As indicated in the previous section it is straightforward to construct the local phase space objects governing 
"reaction" in the NF coordinates, (<7,p). In this section, we will define the various structures in the normal 
form coordinates and discuss briefiy the consequences for the original dynamical system. 

The structure of an energy surface near a saddle point: For < 0, the energy surface consists of 
two disjoint components. The two components correspond to "reactants" and "products." The top panel 
of Fig. [1] shows how the two components project to the various planes of the normal form coordinates. The 
projection to the plane of the saddle coordinates (<7i,pi) is bounded away from the origin by the two branches 
of the hyperbola, qipi = / < 0, where / is given implicitly by the energy equation with the centre actions 
J/c, k = 2, . . . , n, set equal to zero: J^nfI^, 0, . . . , 0) = < 0. The projections to the planes of the centre 
coordinates, (qk^Pk), = 2, . . . , n, are unbounded. 

^The original coordinates (Q, P) typically have an interpretation as configuration space coordinates and momentum coordi- 
nates. The normal form coordinates (q,p), in general, do not have such a physical interpretation since both q and p are nonlinear 
functions of both Q and P. 
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Fi gure ll Projection of energy surfaces (turquoise regions) to the planes of the normal form coordinates. The energy surface 
have energy h < (top panel), h = (middle panel), h > (bottom panel) . 



7 





Figure 2: Projection of the dividing surface and reacting and nonreacting trajectories to the planes of the normal form 
coordinates. In the plane of the saddle coordinates, the projection of the dividing surface is the dark red diagonal line segment, 
which has qi = pi. In the planes of the centre coordinates, the projections of the dividing surface are the dark red discs. 
Forward and backward reactive trajectories (yellow and blue) project to the first and third quadrant in the plane of the saddle 
coordinates, respectively, and pass through the dividing surface. The red and green curves mark nonreactive trajectories on the 
reactant side (pi — qi > 0), and on the product side (pi — qi < 0), of the dividing surface, respectively. The turquoise regions 
indicate the projections of the energy surface. 



At = 0, the formerly disconnected components merge (the energy surface bifurcates), and for > 0, 
the energy surface has locahy the structure of a spherical cylinder, S^^~^ x R. Its projection to the plane 
of the saddle coordinates now includes the origin. In the first and third quadrants it is bounded by the two 
branches of the hyperbola, qipi = / > 0, where / is again given implicitly by the energy equation with 
all centre actions equal to zero, but now with an energy greater than 0: T^TneI^, 0, . . . , 0) = > 0. The 
projections to the planes of the centre coordinates are again unbounded. This is illustrated in the bottom 
panel of Fig. [l] 

The dividing surface, and reacting and nonreacting trajectories: On an energy surface with > 0, 
we define the dividing surface by qi = pi. This gives a (2n — 2)-sphere which we denote by S^g~^(/i). Its 
projection to the saddle coordinates simply gives a line segment through the origin which joins the boundaries 
of the projection of the energy surface, as shown in Fig. [2] The projections of the dividing surface to the 
planes of the centre coordinates are bounded by circles (p^ + g^)/2 = J/e, /c = 2, . . . , where Jk is determined 
by the energy equation with the other centre actions, J^, / 7^ /c, and the saddle integral, /, set equal to zero. 
The dividing surface divides the energy surface into two halves, Pi — ^1 > and Pi — <7i < 0, corresponding 
to reactants and products. 

As mentioned above, trajectories project to hyperbolae in the plane of the saddle coordinates, and to 
circles in the planes of the centre coordinates. The sign of / determines whether a trajectory is nonreacting 
or reacting, see Fig. [2j Trajectories which have / < are nonreactive and for one branch of the hyperbola 
QiPi = I they stay on the reactants side and for the other branch they stay on the products side; trajectories 
with / > are reactive, and for one branch of the hyperbola qipi = / they react in the forward direction, 
i.e., from reactants to products, and for the other branch they react in the backward direction, i.e., from 
products to reactants. The projections of reactive trajectories to the planes of the centre coordinates are 
always contained in the projections of the dividing surface. In this, and other ways, the geometry of the 
reaction is highly constrained. There is no analogous restriction on the projections of nonreactive trajectories 
to the centre coordinates. 

The normally hyperbolic invariant manifold (NHIM) and its relation to the 'activated complex': 

On an energy surface with > 0, the NHIM is given by qi = Pi = 0. The NHIM has the structure of a 
(2n — 3)-sphere, which we denote by S^^m(^)- NHIM is the equator of the dividing surface; it divides 
it into two "hemispheres": the forward dividing surface^ which has qi — pi > 0, and the backward dividing 
surface^ which has qi = Pi < 0. The forward and backward dividing surfaces have the structure of (2n — 2)- 
dimensional balls, which we denote by B'^~^'^{h) and B'^^{h)^ respectively. All forward reactive trajectories 
cross ^dsli^(^)' backward reactive trajectories cross B'^^(h). Since = pi = in the equations of 
motion (Is]) implies that qi = pi = 0, the NHIM is an invariant manifold, i.e., trajectories started in the 
NHIM stay in the NHIM for all time. The system resulting from (71 = pi = is an invariant subsystem with 
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Figure 3: The projection of the NHIM and the local parts of its stable and unstable manifolds, W^{h) and W^{h), to the 
planes of the normal form coordinates. In the plane of the saddle coordinates, the projection of the NHIM is the origin marked 
by the blue bold point, and the projection of W^(h) and W'^(h) are the pi-axis and gi-axis, respectively. W^{h) consists of the 
forward and backward branches Wf{h) and W^(h), which have pi > and pi < 0, respectively; W^{h) consists of W^{h) and 
W^{h), which have gi > and qi < 0, respectively. In the plane of the centre coordinates, the projections of the NHIM, W^{h), 
and W^{h) (the blue circular discs) coincide with the projection of the dividing surface in Fig. |2] The turquoise regions mark 
the projections of the energy surface. 



one degree of freedom less than the full system. In fact, qi = Pi = defines the centre manifold associated 
with the saddle-centre-- • • -centre equilibrium point, and the NHIM at an energy h greater than the energy of 
the quilibrium point is given by the intersection of the centre manifold with the energy surface of this energy 
h [Uzer et al.(2QQ2)| [Waalkens fc Wiggins(2QQ4)| . 

This invariant subsystem with one degree of freedom less than the full system is the "activated complex" 
(in phase space), located between reactants and products. The NHIM can be considered to be the energy 
surface of the activated complex. In particular, all trajectories in the NHIM have / = 0. 

The equations of motion (|8| also show that pi — qi < on the forward dividing surface B'^~^'^(h)^ and 
Pi—qi > on the backward dividing surface B'^^(h). Hence, except for the NHIM, which is is an invariant 
manifold, the dividing surface is everywhere transverse to the Hamiltonian fiow. This means that a trajectory, 
after having crossed the forward or backward dividing surface, B'^~^'^{h) or B'^^(h)^ respectively, must leave 
the neighbourhood of the dividing surface before it can possibly cross it again. Indeed, such a trajectory 
must leave the local region in which the normal form is valid before it can possibly cross the dividing surface 
again. 

The NHIM has a special structure: due to the conservation of the centre actions, it is filled, or foliated^ 
by invariant (n — l)-dimensional tori, T^~^. More precisely, for (i = 3 degrees of freedom, each value of J2 
implicitly defines a value of J3 by the energy equation i^NrCO, ^2, ^3) = h. For three degrees of freedom, the 
NHIM is thus foliated by a one-parameter family of invariant 2-tori. The end points of the parameterization 
interval correspond to J2 = (implying q2 = p2 = 0) and J3 = (implying qs = Ps = 0), respectively. At 
the end points, the 2-tori thus degenerate to periodic orbits, the so-called Lyapunov periodic orbits. 

The stable and unstable manifolds of the NHIM forming the phase space conduits for reactions: 

Since the NHIM is of saddle stability type, it has stable and unstable manifolds, W^{h) and W^{h). The 
stable and unstable manifolds have the structure of spherical cylinders, S^"^~^ x R. Each of them consists of 
two branches: the "forward branches", which we denote by Wj{h) and W^{h)^ and the "backward branches", 
which we denote by W^{h) and W^{h). In terms of the normal form coordinates, Wl{h) is given by = 
with pi > 0, W^{h) is given by pi =0 with (71 > 0, W^{h) is given by gi = with pi < 0, and W^{h) is 
given by pi = with < 0, see Fig. [3] Trajectories on these manifolds have / = 0. 

Since the stable and unstable manifolds of the NHIM are of one less dimension than the energy surface, 
they enclose volumes of the energy surface. We call the union of the forward branches, Wf{h) and W^{h), 
the forward reactive spherical cylinder and denote it by W^{h). Similarly, we define the backward reactive 
spherical cylinder^ Wh{h)^ as the union of the backward branches, W^{h) and W^{h). 

The reactive volumes enclosed by Wf{h) and Wb{h) are shown in Fig. [4] as their projections to the normal 
form coordinate planes. In the plane of the saddle coordinates, the reactive volume enclosed by Wf{h) projects 
to the first quadrant. This projection is bounded by the corresponding hyperbola qipi = /, with / obtained 
from i^NF(^5 • • • 5 0) = h. Likewise, Wi){h) projects to the third quadrant in the (g'i,pi) -plane. Wf{h) 
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Figure 4: Projections of the reactive volumes enclosed by the forward and backward reactive spherical cylinders, W{(h) and 
Wiy(h), and the forward and backward reactions paths, to the planes of the normal form coordinates. The volumes enclosed by 
W{(h) and VKb(/i) project to the dark pink and green regions in the first and third quadrant in the plane of the saddle coordinates, 
respectively. These volumes project to the dark green/dark pink brindled disks in the planes of the centre coordinates, where 
their projections coincide with the projection of the NHIM and the dividing surface in Figs. [2] and [3] The forward and backward 
reaction paths project to the two branches of a hyperbola marked blue in the first and third quadrant in the plane of the saddle 
coordinates, respectively, and to the origins (bold blue points) in the planes of the centre coordinates. The turquoise regions 
mark the projections of the energy surface. 

encloses all forward reactive trajectories; Wi){h) encloses all backward reactive trajectories. All nonreactive 
trajectories are contained in the complement. 

Forward and backward reaction paths: The local geometry of Wf{h) and W\){h) suggests a natural 
definition of dynamical forward and backward reaction paths as the unique paths in phase space obtained by 
putting all of the energy of a reacting trajectory into the reacting mode, i.e., setting ^2 = • • • = = P2 = • • • = 
Pn = 0. This gives the two branches of the hyperbola qipi = /, with / obtained from i^NF(^7 0, . . . , 0) = /i, 
which in phase space are contained in the plane of the saddle coordinates, see Fig. [4j This way, the forward 
(respectively, backward) reaction path can be thought of as the "centre curve" of the relevant volume enclosed 
by the forward (resp., backward) reactive spherical cylinder Wf{h) (resp., W\){h)). These reaction paths are 
the special reactive trajectories which intersect the dividing surface at the "poles" (in the sense of North and 
South poles, where qi = pi assumes its maximum and minimum value on the dividing surface). 

The Transmission Time through the Transition State Region: The normal form coordinates provide 
a way of computing the time for all trajectories to cross the transition region. We illustrate this with a forward 
reacting trajectory (a similar argument and calculation can be applied to backward reacting trajectories). 
We choose the boundary for the entrance to the reaction region to he pi — qi = c for some constant c > 0, 
i.e., initial conditions which lie on the reactant side of the transition state, and the boundary for exiting the 
reaction region to he pi — qi = —c on the product side. We now compute the time of flight for a forward 
reacting trajectory with initial condition on pi — qi = c to reach pi — qi = —c on the product side. The 
solutions are qi{t) = qi{0) exp(A(/, J2, . . . , Jn)t) and pi{t) = pi{0) exp(— A(/, J2, . . . , Jn)t) (see ([8|), where 
A(/, J2, . . . , Jn) is determined by the initial conditions. This gives the time of fiight as 

r=(A(/,J2,...,J„))-^ln(^^). (11) 

The time diverges logarithmically as qi{0) 0, i.e., the closer the trajectory starts to the boundary Wf{h). 
It is not difficult to see that the time of ffight is shortest for the centre curve of the volume enclosed by Wf(/i), 
i.e., the trajectory which traverses the transition state region fastest is precisely our forward reaction path. A 
similar construction applies to backward reactive trajectories. 

In fact, the normal form can be used to map trajectories through the transition state region, i.e. the 
phase space point at which a trajectory enters the transition state region can be mapped analytically to the 
phase space point at which the trajectories exits the transition state region. 
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2.3 The foliation of the reaction region by Lagrangian submanifolds 



The existence of the n integrals of motion, (/, J2, . . . , Jn), induce phase space structures which lead to further 
constraints on the trajectories in addition to the ones described above. In order to describe these structures 
and the resulting constraints it is useful to introduce the so called momentum map, M |Guillemin(1994)| 
Marsden fc Ratiu(1999)| which maps a point (<7i, . . . , ^n^Pi, • • • ^Pn) in the phase space x to the n 
integrals evaluated at this point: 

: X ^ , (gi, . . . , . . . ^ (/, J2, . . . , Jn) . (12) 

The preimage of a value for the constants of motion (/, J2, . . . , Jn) under Ai is called a fibre. A fibre thus 
corresponds to the common level set of the integrals in phase space. 

A point (g'l, . . . , g'n^Pi, • • • ^Pn) is called a regular point of the momentum map if the linearisation of the 
momentum map, DM, has full rank n at this point, i.e., if the gradients of the n integrals /, J/e, /c = 2, . . . , n, 
with respect to the phase space coordinates {q,p) are linearly independent at this point. If the rank of DAi 
is less than n then the point is called irregular. A fibre is called regular if it consists of regular points only. 
Else it is called an irregular fibre. In fact almost all fibres are regular. They are n-dimensional manifolds 
given by the Cartesian product of an hyperbola qipi = / in the saddle plane {qi,Pi) and n — 1 circles in 
the centre planes A: = 2, . . . , n. Since the hyperbola qipi = I consists of two branches each of which 

have the topology of a line R, the regular fibres consist of two disjoint toroidal cylinders, T"^~^ x R, which 
are the Cartesian products of a (n — l)-dimensional torus and a line. We denote these toroidal cylinders by 

A|,j„...,j„ = {iq,p) e ]R2" : pm =I,l{pl+ql)=J2,-..,l {pI + ql) = J„ , 9i > 0} (13) 

and 

A7,J„...,J. ={fep) ^I^'" • PlQl=I. ^{pI^qD =^2,... ,-(Pn+^n) =^n,^l <0}. (14) 

A^j^ Jd J2 Jd Lagrangian manifolds [Arnold(1978)| . The Lagrangian manifolds consists of all 

trajectories which have the same constants of motion. In particular the Lagrangian manifolds are invariant, 
i.e. a trajectory with initial condition on a Lagrangian manifold will stay in the Lagrangian manifold for all 
time. For / < 0, the Lagrangian manifolds j and ^'j j consist of nonreactive trajectories in the 

react ants resp. products components of the energy surface. For / > 0, consists of forward reactive 

trajectories, and Aj consists of backward reactive trajectories. 

We will discuss the irregular fibres in more detail in Sec. [4j 



2.4 Imphcations for the original system 

The normalization procedure proceeds via formal power series manipulations whose input is a Taylor expan- 
sion of the original Hamiltonian, H, necessarily up to some finite order, M, in homogeneous polynomials. 
For a particular application, this procedure naturally necessitates a suitable choice of the order, M, for the 
normalization, after which one must make a restriction to some local region, £, about the equilibrium point in 
which the resulting computations achieve some desired accuracy. Hence, the accuracy of the normal form as a 
power series expansion truncated at order M in a neighborhood C is determined by comparing the dynamics 
associated with the normal form with the dynamics of the original system. There are several independent 
tests that can be carried out to verify accuracy of the normal form. Straightforward tests that we use are 
the following: 

• Examine how the integrals associated with the normal form change on trajectories of the full Hamilto- 
nian (the integrals will be constant on trajectories of the normal form). 

• Check invariance of the different invariant manifolds (i.e. the NHIM and its stable and unstable 
manifolds) with respect to trajectories of the full Hamiltonian. 

Both of these tests will require us to use the transformations between the original coordinates and the 
normal form coordinates. Specific examples where M, C and accuracy of the normal forms are considered can 
be found in [Waalkens et al.(2QQ4a)[|Waalkens et al.(2QQ4b)'||Waalkens et al.(2QQ5a)[|Waalkens et al.(2QQ5b)[ 



Waalkens et al.(2QQ5c)[ T 
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3 Models of the local phase space structures for 2-DoF systems 



In the discussion of the phase space structures governing reaction type dynamics in Sec. \2.2\ we visuahzed 
these structures by projecting them to the 2-dimensional planes of the coordinates and conjugate momenta 
of the normal form. The topology of the phase space structures and the interrelationship between them is 
to a large extent obscured in these projection because of the lacking injectivity of the projections. For n = 2 
DoFhowever, a surface of constant energy containing the phase space structures is three-dimensional. For 
this case, it is therefore possible to visualize the phase space structures in three-dimensional space. In this 
section we will develop three models to explicitly implement such visualizations. Our discussion will be in 
the context of linear Hamiltonian systems since they allow us to describe the phase space structures in terms 



of explicit formulae. The influence of nonlinear terms is discussed in Sec. 3.4 
Therefore we begin by considering the quadratic Hamiltonian 

with A, > 0. 

This gives the linear Hamiltonian vector field 



(15) 



X = 

Px 

y = 

Py = 



dH 

dpx 
dH 

dx 
dH 

dpy 

dH 
dy 





(16) 


Xx, 


(17) 


OjPy, 


(18) 


-uy, . 


(19) 



The Hamiltonian (|T5|) is already in normal form. To keep the notation simple we will not specifically 

(20) 



indicate this by a subscript NF. In fact we have 

H{x, y.px.py) = K{I, J) = XI ^ LuJ 



where / and J are the integrals given in (|6|) and ([7|), respectively. 

For a fixed energy /i, the three dimensional energy surface is given by the level set of the Hamiltonian (15) 
as follows: 



^{h) = {{x,y,p^,Py) gR^ : X{pI-x^) 
We rewrite the energy equation in the form 

Xpl + Ljpl + Ljy'^ = 2/i + Xx^ , 



2h}, 



(21) 



(22) 



and fix X such that the right hand side of this equation is positive. Note that if h > the right hand side 
will be positive for any fixed x G R, while if < it will be positive only for \x\ > ^/—2h/X^ i.e. for x 
sufficiently small or large. For such an Eq. (22) defines an ellipsoid in the three-dimensional Euclidean 
space, R^, with coordinates (Px^Py^y)- This ellipsoid is rotationally symmetric about the p^^-axis. It has two 

1/2 

semi-axes of length ((2/i + Xx'^)/u) and one semi-axis of length ((2/i 
a two-dimensional sphere, S^, which we denote by (/?'), i.e. 



1/2 

Ax^)/A) . Topologically, this is 



E,l{h) = {{p,,Py,y)e-R^ : Xpl + u;{pl + y^) = 2h + \x^} . 



(23) 



The energy surface thus appears to have the structure of a one-parameter family of two-dimensional spheres 



parametrized by x, where the 'radius' (the length of the shortest semi- axis of the ellipsoid (22)) can shrink to 
zero if h < 0. The three models of the energy surface and the phase space structures contained in it that we 
present in the following differ by the way of representing the 2-spheres S^(/i). For clarity, let us summarize 
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the phase space structures discussed in Sec. [5] for the Hamiltonian (15) before we proceed. We wih restrict 
ourselves to the case h > 0. 

The dividing surface is obtained by setting x = on the energy surface. This gives the two-dimensional 
sphere 



^l,{h) = n{x = 0} = {{x, y,p,,py) e : Xpl + + y^) = 2h, x = 0} . 



(24) 



For the 2-DoF case, the NHIM is an unstable periodic orbit (the Lyapunov periodic orbit associated with the 
saddle point). For the Hamiltonian (15), it is obtained by setting x = and = on the energy surface. 
This gives the circle 

Snhim(^) = n {x = 0, = 0} = {(x, y,p:,,py) e : u;{pl + y^) = 2h, x = 0, = 0} . (25) 
The NHIM divides the dividing surface, S^g(/i), into the forward hemisphere 



Blj{h) = ^{h) n {x = 0, > 0} = {{x,y,p,,py) G : = ^(2/i -u{pl + y^))/X, x = 0} (26) 
and the backward hemisphere 



Bl^^{h) = S(/i) n {x = 0, < 0} = {(x, y,P:,,Py) G : = {2h - u;{pl ^ y^)) /X, x = 0} , (27) 

which are two dimensional balls or disks. The stable manifold of Sj^jjjj^(/i) consists of the forward branch 

Wf = ^{h) n{x = x<0} = {{x,y,p^,py) e : co{pl + y^) = 2h, x = x < 0} (28) 
and the backward branch 

= ^{h) n{x = -p^, x>0} = {(x, y,P:,,Py) e : u;{pl + y^) = 2h, x = -p^, x > 0} . (29) 
Similarly, the unstable manifold of Snhim(^) consists of the forward branch 

W^^ = S(/i) n{x=p,,x>0} = {{x,y,p,,py) e R^ : + ^') = 2h, x = p,, x > 0} (30) 

and the backward branch 

= S(/i) n{x=p,,x<0} = {(x, y,p^,py) eR^ : u{pI^ y^) = 2h, x = p„ x < 0} . (31) 

Finally, the forward and backward reaction paths are obtained from putting all the energy in the x DoF. 
This gives the lines 



S(/i) n{y = py=0,p:,>0} = {{x,y,p^,py) G R^ : = V(2/i + Ax)/A, y = Py = 0} (32) 
for the forward reaction path, and 

S(/i) n{y=py = 0,p:,<0} = {(x, y,P:,,Py) G R^ : = - V(2/i + Ax)/A, y = py = 0} (33) 
for the backward reaction path. 



Note that the energy surface (|21|, the phase space structures (|24[)-(|33[) and also the Lagrangian cylinders 

(34) 



defined in (13) and (14) are invariant under the symmetry action 

{x,y,Px,Py) {x,e~''^y,p^,e~''^py) , (p eE>^ 
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3.1 Model I: projection of 2-spheres to their equatorial planes 



The main idea of the first model is to view the 2-sphere (23) as consisting of two hemispheres and an equator, 
i.e. as the disjoint union of two open two-dimensional disks or balls, 5^, and a circle, or equivalently one- 
dimensional sphere, S^, and using parametrizations of these geometric structures to parametrize the energy 
surface. 

The decomposition of a 2-sphere, S^, in terms of 2-balls and a circle is familiar from visualizing the surface 
of the Earth by drawing two filled-in disks (5^), one each for the Northern and Southern hemispheres, and 
identifying the boundaries of the two disks, 

= ^North U ^Equator ^ ^South^ (35) 

which we illustrate in Fig. |5] 




Figure 5: A model of the sphere formed by taking two open discs (2-balls) ^North ^South 
identifying points on their boundaries; the "equator" S^quator- ^^^^ shown is the representation of a smooth 
closed curve on which passes through the North and South poles, N and S; the points at which the 
representation of the curve "jumps" from one disc in the model to the other disc are marked by dots. 



For the 2-spheres S^(/i) in Eq. ([23|), we define the decomposition (35) as 



Sl{h) = Bl,{h)uSlJh)UBl^,{h), 



(36) 



where 



Bljh) = {{p:c,Py,y) e : co{pI + < 2/1 - Xx^, = ^{2h + Xx^ - upl - ujy^)/X > 0} , (37) 
B-Jh) = {{P:c,Py,y) e : co{pI + ?/2) < 2/1 - Xx^ p, = -^{2h + Xx^ - wp2 _ u;y^)/X < 0} , (38) 
Sj.^W = {{Px,Py,y)eR^ ■■ iv{pl + y'')=2h-Xx\p^ = 0}. (39) 



We can view the manifolds B\ ,j.{h) and i?i^(/i) as graphs of the functions 



P±,, : B2(/i) ^ R , {y,py) ^ ±^{2h + Ax^ - Lopl - Loy^)/X , 
where the domains are the circular disks 



(40) 



Bl{h) = {{py,y)eR^ : u;{pl^y^) <2h^Xx^}cR\ (41) 

The graphs 'join' at the equator xW- order to geometrically visualize S^(/i) we take two copies of the 
circular disks B'^{h) in the (y,^^/) -plane where Px is given by the function Vj^^x on one copy which we denote 
by B\ .j.{h) and the function V-^x on the other copy which we denote by B'^ .j.{h). B'^ .j.{h) thus denote the 
images of B'j^ .j.{h) under the projection 



^P.x • B^^^{h) , {Px^Py^y) ^ {y^Py) , 



(42) 
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Figure 6: Model I: projection of S^(/i) to the (?/,p^)-plane. 

see Fig. [oj Along the boundaries of B'^^^{h) and B'^^^{h) we have Px = 0, and the points on the boundaries 
which have the same {y^Py) are identified. 

The first model (Model I) of the energy surface and the phase space structures contained in it is obtained 
from defining a map from the energy surface ^{h) to using the projections TTp^x^ x G R, in (42) to map 
points {x^y^Px^Py) G ^{h) according to 

: {x,y,Px,Py) ^ {x.TTp^xiPx^Py^y)) = {x,y,py). (43) 



The map is singular along the cylinder given by the union of the circles Sj ^{h) in (39) with x G R. This 
singular cylinder is mapped to the quadric 

■ r (44) 



2h/u; 2h/X 2h/uo 

which is a single sheeted hyperboloid that is rotationally symmetric about the x axis. Each point in 
(x, y^Py) G with 

' <1 (45) 



2h/uj 2h/\ 2h/u 

has two preimages that differ by the sign of Px- We can thus represent the energy surface T,{h) in the three 
dimensional space, R^, in terms of the manifolds 



z/2 (46) 
where Px = (2/i + x^ ~Py~ iP'Y^^ positive on Y}^(K) and Px = —(2h^x^ —py — y^Y^^ negative on Y}_(K). 



The full energy surface S(/i) is obtained by gluing S^^(/i) along their boundaries (44) where Px = 0. We show 
illustrations of S^(/i) and Ti^_{h) in the top panels of Fig. ^ As a consequence of the symmetry ( [34| ) the 
components ^±{h) are both rotationally symmetric about the x-axis. It is therefore useful to consider also 
the intersections of ^\{h) and ^\{h) (and of the other phase space structures discussed in the following) with 
the plane = as shown in the bottom panels of Fig. [7| It is easy to see the wide-narrow- wide geometry of 
the bottleneck type structure of the energy surface in this representation of the energy surface. The reactants 
region is on the left, the product regions is on the right. 



^We note that for /i < the quadrics (|44|> are two-sheeted hyperboloids. For h = 0, the quadric is a 'diabolo'. 
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Figure 7: Model I: the energy surface. The orientation of the x, y and Py axes is indicated by the dreibein. 
(Note that the origin of the coordinate system is not the center of the dreibein.) The top left panel shows 
Tj\(h)] the top left panel shows T}_{h). The two panels at the bottom show the intersections of Tj\{h) and 
SL(/?') with the plane y = 0. 



16 



Figure 8: Model I: nonreactive trajectories that stay in the reactants component. The red wireframes in the 
top panels indicate the two pieces of a Lagrangian cylinder Aj j defined in (14) for /< contained in ^\{h) 
and SL(/i^), respectively. The black line shows one of the trajectories contained in Aj j. The red lines in the 
bottom panels show the intersection of Ajj with the plane y = 0. 
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Figure 9: Model I: nonreactive trajectories that stay in the products component. The green wireframes in 
the top panels indicate the two pieces of a Lagrangian cylinder A.'\ j defined in (13) for /< contained in 
T}j^(h) and Y}_(h)^ respectively. The black line shows one of the trajectories contained in A]^j. The green 
lines in the bottom panels show the intersection of K'\ j with the plane ^ = 0. 



In Fig. [s] we show the image under of a Lagrangian cylinder Aj j defined in (14) for /< together 
with a nonreactive trajectory contained in this cylinder. Note that for a 2-DoF system and given values for 
the integral / and the energy /i, the second integral J is fixed by the energy equation J) = h (see (20)). 
The trajectory approaches the bottleneck region with positive Px from the react ants component of the energy 
surface within ^\{h) (the left panels of Fig. p|, reaches the boundary (44) of S+ (/?-), where Px = and the 



trajectory 'jumps' to ^L{h) (the right panels in Fig.|8| after which it returns back deeper into the reactants 
region withpa^ < 0. 

Figure 9 shows the image of a Lagrangian cylinder Afj for / < together with a trajectory contained 
in A'^ J. The trajectory approaches the bottleneck region with negative Px from the products side of the 
energy surface within SL(/i) (the right panels of Fig. [9|, reaches the boundary (44) of (/?-), where Px = 
and the trajectory 'jumps' to ^\{h) (the left panels in Fig. |9|, in which it returns back deeper in to the 
products region with Px > 0. 



Figure 
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shows the F^ images of Lagrangian cylinders Aj j and A'f j for / > 0. For / > 0, A'f j contains 
forward reactive trajectories which have Px > and therefore are completely contained in ^\{h) (the left 
panels in Fig. 10), and Ajj contains backward reactive trajectories which have Px < and therefore are 
completely contained in ^^_{h) (the right panels in Fig. 10). The forward and backward trajectories spiral 



about the forward and backward reaction paths which are also shown in Fig. 10 
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The images of the dividing surface, S^g(/i), and the NHIM, Snhim(^)' under F^ are shown in Fig. 
The forward hemisphere B'^^ ^{h) is contained in ^\{h) (the left panels in Fig. 11); the backward hemisphere 
^ds b(^) contained in ^^_{h) (the right panels in Fig. 11). The periodic orbit Sj^HiM(^) which forms the 



18 




Figure 10: Model I: forward and backward reactive trajectories spiraling about the forward and backward 
reaction paths, respectively. The yellow wireframe in the top left panel indicates a Lagrangian cylinder K'\ j 
with / > contained in T}^{h)] the blue wireframe in the top right panel indicates a Lagrangian cylinder 
A^j with / > contained in T}_(h). The black lines show trajectories contained in the respective cylinders. 
The blue centerlines show the forward (left) and backward (right) reaction paths. The bottom panels show 
the intersections of the Lagrangian cylinders with the plane ^ = together with the forward and backward 
reaction paths. 
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Figure 11: Model I: the dividing surface and the NHIM. The dark red wireframes indicate the two hemispheres 
^ds, f(^) (^^P ^^^^ panel) and ^^s, b(^) (^^P right panel) the dividing surface E>'^^{h). The light blue circle in the 
top panels shows the periodic orbit Sj^HiM(^)- lower panels show the intersections of ^(/i), B'^^ ^{h) 
and Sjjjjj|^(/i) with the plane y = 0. 
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Figure 12: Model I: the stable and unstable manifolds of the NHIM. The wireframes in the top left panel 
show the branches Wf (green) and W^^ (red), and the wireframes in the top right panel show the branches 
(green) and (red) of the stable and unstable manifolds of Snhim(^)- edich branch, a trajectory 
(black line) contained in the respective branch is shown. The bottom panels show the intersections of the 
branches W^/^^ with the plane y = 0. 



equator of SL(/i) has Px = and thus runs along the boundaries of ^\{h) and SL (/?-). 

Figure [l2| shows the images of the stable and unstable manifolds and W"^ of Sjjfjjj^(/i). The forward 

), and the backward branches and 
lows how trajectories contained in the 



branches Wf and W}^ are contained in 5]^(/i) (left panels in Fig. 



are contained in T^L{h) (left panels in Fig. 12). Figure 



branches W^!^ asymptotically approach the periodic orbit §nhim(^) respective time limit t ±oo. 
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also s 



^f/b 



Regarding Fig. 12 together with Fig. 10 we see that the forward reactive cylinder Wf U W}^ contained in 
^\{h) encloses all forward reactive trajectories, and the backward reactive cylinder UW{^ encloses all the 
backward reactive trajectories. Note that the images of the forward and backward reactive cylinders are 
smooth along their junctions given by the periodic orbit SJ^him(^) while in the original phase space WfUW^ 
and U have a kink along SJ^himI^)- This is an artifact of the map F which, as mentioned above, is 
singular along the cylinder U^^^iaSo ^{h) which contains Snhim(^) ~ o(^)- 



3.2 Model II: stereographic projection of 2-spheres 



The main idea of the second model (Model II) is to represent the 2-spheres S^(/i) defined by Eq. (23) in terms 
of stereographic type projections to a plane. Recall that S^(/i) is given by the axisymmetric ellipsoid (23). 
We define a stereographic projection of this ellipsoid by projecting each point on it along the line through 
this point and the South pole (px^Py^y) = (— \/(2/i + Ax^)/A, 0, 0) of the ellipsoid to the equatorial plane 
of the ellipsoid. This projection is illustrated in Fig. 13 Points on the hemisphere ^{h) of (/?-), which 



we defined in (37), are projected to the circular dislTbounded by the equator ([39]) (e.g., the point Pi in 
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Figure 13: Model II: stereographic projection of S^(/i) to its equatorial plane with coordinates {y^Py). The 



points Pi and P2 are the images of the points Pi and P2 under the projection (47). 



Fig. 13). The North pole is mapped to the center of this disk. Points on the hemisphere B_ ^{h) defined in 



(38) are projected to points in the plane outside of the disk (see, e.g., the point P2 in Fig. [13]). The South 
pole {pxiPyiV) = ( — \/ (2/1 + Ax^)/A, 0, 0) is mapped to infinity. The equator ^^^^(h) of S^(/?') (see ([39|) is 
invariant under the projection. Formally this projection is given by 



{Px^Py^y) ^ {y^Py) : 



where 



y^{2h + Xx^)/\ 
v/(2/i + Aa;2)/A+p 



Pv 



Py^{2h + \x^)/\ 
v/(2/i + Aa;2)/A+p^ 



(47) 



(48) 



and R denotes the compactified plane R? U {00}. Note that the familiar standard stereographic projection 
is from a 2-sphere of unit radius to its equatorial plane where the 'base point' of the projection is chosen to 
be the North pole rather than the South pole like in our choice. This explains the scaling factors and signs 
in (|48| which differ from the standard formulae. We use the projection tt^^^ in (47) to define a map of points 
{x^y^Px^Py) on the energy surface T,{h) according to 

P" : {x,y,Px,Py) ^ {x.TTsgxiPx^Py^y)) = {x,y,py). (49) 

The map P^^ is singular along the North poles of the spheres (/?-), x G R, where y and Py are zero, and 

Px reaches its minimal value for fixed x and h. The set of singular points thus is the backward reaction 

2 

path which is mapped to infinity. As a result the image of Ti(h) under P^^ is the Cartesian product R x R . 
Note that the map P^^ preserves the symmetry (34), or more precisely, like in Model I, manifolds which are 



invariant under the symmetry action (34) are rotationally symmetric about the x axis of the image of P^^. 

1 visualizations of the various phase space 

2 

Since the energy surface maps to R x R 



In the following we use this representation to show geometrical visualizations of the various phase space 

I I ^ _2 

structures that govern reaction dynamics analogously to Sec. 



3.1 



under P^^ we omit a picture of the image of the energy surface. In Fig. 14 we show the P^^ image of a 
Lagrangian manifold Aj j with / < together with a nonreactive trajectory contained in this cylinder. The 
trajectory comes in from the reactants component (x <C —1) on the left with positive Px^ spiraling about the 
X axis. The radius of the rotation increases as x increases. After reaching a maximal x value (where Px = 0) 
the trajectory returns deeper and deeper into the reactants region, still spiraling about the x axis with the 
radius of rotation increasing as x decreases. In Fig. 15 we show the analogous picture for the P^^ image 



of a Lagrangian manifold A'^ j with / < 0. The trajectories come in from the products component on the 
right (x > 1) with negative Px undergoing rotations with increasing radius about the X axis as x decreases. 
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Figure 14: Model II: nonreactive trajectories that stay in the reactants component. The orientation of the 
X, y and Py axes is indicated by the dreibein. The wireframe in the left panel shows a Lagrangian cylinder 
A J J defined in (14) for / < 0. The black line shows one of the trajectories contained in Aj j. The red lines 
in the right panel show the intersection of Ajj with the plane ^ = 0. 




Figure 15: Model II: nonreactive trajectories that stay in the products component. The green wireframe in 



the left panel shows a Lagrangian cylinder A'f j defined in (13) for / < 0. The black line shows one of the 
trajectories contained in A'f j. The green lines in the right panel show the intersection of A'^ j with the plane 
^ = 0. ' 
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Figure 16: Model II: forward reactive trajectories spiraling about the forward reaction path. The yehow 
wireframe in the left panel indicates a Lagrangian cylinder K'\ j with / > 0. The black line shows a trajectory 
contained in A]^j. The blue line shows the forward reaction path. The right panel shows the intersection of 
J with the plane ^ = together with the forward reaction path. 




Figure 17: Model II: backward reactive trajectories. The blue wireframe in the left panel indicates a 
Lagrangian cylinder Aj j with / > 0. The black line shows a trajectory contained in Aj j. The right panel 
shows the intersection of A^ j with the plane ^ = 0. 



After reaching a minimal x value, the trajectory returns deeper and deeper into the products region, still 
undergoin g ro tation about the x axis, where the radius of the rotation now increases as x increases. 

shows the F^^ image of a Lagrangian cylinder Aj j with / > together with the forward 



Figure 



16 



reaction path which coincides with the x axis. Trajectories contained in K'\ j evolve from the react ants to 
the products spiraling about the forward reaction path. Figure 17 shows the F^^ image of a Lagrangian 
cylinder Aj j with / > that contains backward reactive trajectories. These trajectories spiral with much 
larger radii of rotation than the forward reactive trajectories in Fig. 16 In fact, since the backward reaction 
path is mapped to the line RU oo the backward reactive trajectories can again be considered to spiral about 
the backward reaction path. 

In Fig. 18 we show the images of the dividing surface ^l^{h) and the NHIM SJ^himI^) ™<^er F". The 
dividing surface ^\^{h) maps to the compactified plane x = 0. The F^^ image of the NHIM, §nhim(^)' divides 
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Figure 18: Model II: dividing surface and NHIM. The dark red wireframe in the left panel shows the 
dividing surface S^g(/i). The periodic orbit SJ^him(^) (blue circle) divides it into the forward hemisphere 
^ds f(^) (iiiiier circular disk) and the backward hemisphere B'^^ ^{h) (outer infinite annulus). The right 
panel shows the intersection of '^'^^{h) and SJ^him(^) ^^^^ plane y = 0. 




Figure 19: Model II: forward reactive cylinder, forward dividing surface and NHIM. The green and red 
wireframes in the left panel indicate the branches Wf and W}^ of the stable resp. unstable manifolds of the 
NHIM Sjjjjj|^(/i) (blue circle). The black lines show trajectories contained in Wf and W^. The dark red 
wireframe at the center which is difficult to see shows the forward hemisphere ^{h) of the dividing surface. 
The right panel shows the intersections of the structures on the right with the plane ^ = 0. 



it into the forward hemisphere B'^^ ^{h) which appears as a circular disk in the plane x = centered about 
the origin and the backward hemisphere ^^s, b(^) which appears as an annulus that extends to the boundary 



of at infinity. The F^^ image of the forward reactive trajectories in Fig. 16 which are located near the x 
axis cross the F^^ image of B'^^ ^(/i), and the F^^ image of the backward reactive trajectories in Fig. |T7| which 
stay away from the x axis cross the F^^ image of B'^^ b(^)- "^^^ types of trajectories are enclosed by the 
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forward and backward reactive cylinders Wf U and U whose F^^ images are shown in Figs, 
and [20} respectively. 

Note that in this geometrical representation of the energy surface (Model II) the bottleneck area of Ti{h) 
associated with the reaction region containing the dividing surface is not so pronounced as in Model I. 
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Figure 20: Model II: backward reactive cylinder, backward dividing surface and NHIM. The green and red 
wireframes in the left panel indicate the branches and of the stable resp. unstable manifolds of 
the NHIM SJjhim(^) (blue circle). The black lines show trajectories contained in and W^. The dark 
red wireframe shows the backward hemisphere B'^^ ^{h) of the dividing surface. The right panel shows the 
intersections of the structures on the right with the plane y = 0. 



3.3 Model III: McGehee representation 

The main idea of the third model (Model III) of the energy surface S(/i) is to represent the family of 



(topological) 2-spheres, (/?-), in Eq. (23) parametrized by x G IR as a 'nested' set of concentric 2-spheres in 
R^. This can be achieved by projecting each of the 2-spheres, S^(/i) to a 2-sphere of a radius which depends 
on X. One way to realize this is to map a point P = {px^Py, y) on the S^(/i) according to 



TTMcGx ' {y,Px,Py) ^ {y,Px,Py)l 



s{x) 



--' (y^Px^Py) . 



\\{y^Px^Py)\\ 

Here ||-|| denotes the Euclidean norm, i.e. IKp^^P?/, y)|| = {Px + iP'Y^^ ^i^nd s is the scaling function 



(50) 



5(x) := 1 + - tanhx (51) 

which maps the domain x G IR monotonicly to the interval (1/2, 3/2). The projection TTMcGa^ is illustrated in 
Fig. [21] We name this projection after McGehee who first used this type of construction to present an energy 
surface of the structure x H [McGehee(1969)| . A similar construction can also be found in the work by 

to define a map F^^^ of points {x^y^Px^Py) on the energy surface T^{h) to 

F"^ : {x,y,Px,Py) ^ TTMcGxiPx^Py^y) {y.Px^Py)- (52) 

Like the maps and the map F^^^ also preserves the symmetry. In Model III, manifolds that are 
invariant under the symmetry action have F^^^ images that are symmetric about the Px axis. 

The image of the energy surface ^{h) of topology R x under the map F^^^ is the 'spherical shell' 
consisting of the open ball of radius 3/2 in minus the closed ball of radius 1/2. This image is shown in 

1 is contained in the outer region of this spherical shell; the products 
we also marked the points {y^Py^Px) = (0^ 0^^1/2) 



MacKay MacKay(1990)'. 

We use ttmcGx in iSo 
according to 
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Fig. [22] The reactants region x <C 
region x 1 is contained in the inner region. In Fig 

and (y^Py^Px) = (0, 0, ^3/2) to which the F^^^ images of all trajectories converge as time goes to — oo or +oo 
(as long as the trajectory is unbounded in the respective limit). This is a consequence of the conservation of 
the integrals / and J (the energy in the x and y degrees of freedom, respectively) which implies that y and 



Py are bounded for each trajectory whereas Px diverges as x goes to — oo or +oo. In Fig. 23 we show the 
example of the F^^^ image of a Lagrangian cylinder Ajj defined for / < 0. The trajectories in Ajj come in 
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Figure 21: Model III: McGehee projection of a point P = (y^Px^Py) on the (topological) 2-spliere S^(/i) 
(inner dark shaded sphere) to the point P = P = (y^Px^Py) on the sphere of radius s{x) (outer light shaded 
sphere) with s{x) defined in (51). Note that, depending on the outer sphere could also be inside of the 
inner sphere, or the two could overlap. 



from deep in the reactants region where x <C — 1 and Px ^ 1 which is asymptotically (in the limit t ^ - 
mapped to (y^Py^Px) = (0? 0^3/2) under F^^^ and return deep into the reactants region where x <C — 1 
Px ^ ~^ which is asymptotically (in the limit t oo) mapped to {y^Py^Px ) = (0, 0, -3/2) under F^^^ 



-oo) 
and 



In Fig. 
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we show the F^^^ image of a Lagrangian cylinder A^j with / < which contains nonreactive 
trajectories m the products region. They come in from deep in the products region x 1 with Px <^ —I 
which under F^^^ is mapped to (y^Py^Px) = (0, 0, ~l/2) and move back deep into the products region x ^ 1 
with px which under F^^^ is mapped to (y^Py^Px) = (O^O, 1/2)- 



Fig. 

consists 
Fig. 
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25 
Fig. 



shows the F^^^ image of a Lagrangian cylinder A^j which has / > 0. The Lagrangian manifold 
forward reactive trajectories which spiral about the forward reaction path that is also shown in 

which has / > and hence consists 
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shows the analogous image of a Lagrangian cylinder A^ j 
of backward reactive trajectories. 

For clarity we show the F^^^ images of the branches Wf, 



and of the stable and unstable 
The F^^^ image of the union Wf U W^^ forms the forward 
reactive cylinder in Model III. It encloses the^F^^^ image of all forward reactive trajectories as can be seen 

Similarly the F^^^ image of the backward reactive cylinder U 



manifolds of the NHIM separately in Figs. 28pT 



from comparing Fig. 25 with Fig. 32 



encloses the F^^^ image of all backward reactive cylinder as can be seen from comparing Fig. 26 with Fig. 33 



3.4 Implications for Nonlinear Hamiltonian Vector Fields 

We have chosen to develop the three models of the phase space structures near a saddle x centre equilibrium by 
considering a quadratic, two degree of freedom Hamiltonian which is the normal form, through second order, 
in the neighborhood of such an equilibrium point. The advantage of this is that is provided very explicit 
formulae for the phase space structures that allowed straightforward graphical representation. However, 
precisely the same geometrical structures exist for nonlinear Hamiltonian vector fields in the neighborhood of 
a saddle x center equilibrium point. In the nonlinear case the phase space structures can again be described 
by explicit mathematical formulae using a normal expansion at the saddle x center equilibrium (see Sec. |2.2 



and ll Waalkens et al.(2QQ8)| for a detailed discussion). For a given system, this in general leads a Hamilton 
function which is a power series rather than a linear function of the integrals / and J (see ([2Q|)). The 
coefficients of this powers series are obtained from the normal form calculations for an explicit systems. 
Nevertheless, the normal form theory allows us to conclude that the same geometrical structures, having 
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Figure 22: Model III: the energy surface. The left panel shows the energy surface in the threedimensional 
space with coordinates {y^Py^Px 

). The orientation of the y, Py and Px axes are indicated by the dreibein. 
(Note that the center of the coordinate system is not the center of the dreibein.) The black points mark the 
points {y^Py^Px) = (0, 0, ±1/2) and {y^Py^Px) = (0, 0, ±3/2) (see text). The right panel shows the intersection 
of the energy surface with the plane ^ = 0. 




Figure 23: Modell III: nonreactive trajectories that stay in the reactants component. The red wireframe in 
the left panel indicates a Lagrangian cylinder Ajj defined in (14) for / < 0. The black line shows one of the 
trajectories contained in Aj j. The red lines in the right panel show the intersection of Ajj with the plane 
^ = 0. ' 
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Figure 24: Modell III: nonreactive trajectories that stay in the products component. The green wireframe 
in the left panel indicates a Lagrangian cylinder K'\ j defined in (13) for / < 0. The black line shows one of 
the trajectories contained in A]^j. The green lines in the right panel show the intersection of K'\ j with the 
plane ^ = 0. 




Figure 25: Model III: forward reactive trajectories spiraling about the forward reaction path. The yellow 
wireframe in the left panel indicates a Lagrangian cylinder K'\ j with / > 0. The black line shows a trajectory 
contained in A]^j. The blue line shows the forward reaction path. The right panel shows the intersection of 
K'\ J with the plane ^ = together with the forward reaction path. 
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Figure 26: Model III: backward reactive trajectories spiraling about the backward reaction path. The blue 
wireframe in the left panel indicates a Lagrangian cylinder Aj j with / > 0. The black line shows a trajectory 
contained in Aj j spiraling about the backward reaction path marked in blue. The right panel shows the 
intersection of Aj j with the plane ^ = together with the backward reaction path. 




Figure 27: Model III: dividing surface and NHIM. The dark red wireframe in the left panel shows the 
dividing surface S^g(/i). The periodic orbit SJ^him(^) (blue circle) divides it into the forward hemisphere 
^ds f(^) (^^^ hemisphere on the right) and the backward hemisphere B'^^ ^{h) (the hemisphere on the left). 
The right panel shows the intersection of '^'^^{h) and Snhim(^) ^^^^ plane ^ = 0. 
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Figure 28: Model III: the forward branch of the stable manifold. The green wireframe in the left panel 
shows Wf. The black line shows a trajectory contained in Wf . The green lines in the right panels show the 
intersection of Wf with the plane ^ = 0. 




Figure 29: Model III: the forward branch of the unstable manifold. The red wireframe in the left panel 
shows W^. The black line shows a trajectory contained in W^^. The red lines in the right panels show the 
intersection of with the plane ^ = 0. 
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Figure 30: Model III: the backward branch of the stable manifold. The green wireframe in the left panel 
shows W^. The black line shows a trajectory contained in W^. The green lines in the right panels show the 
intersection of with the plane ^ = 0. 




Figure 31: Model III: the backward branch of the unstable manifold. The red wireframe in the left panel 
shows W^. The black line shows a trajectory contained in W^. The red lines in the right panels show the 
intersection of with the plane ^ = 0. 



32 



Figure 32: 
of Figs. [27) 



Model III: the forward reactive cylinder, the dividing surface and the NHIM as a superposition 
281 and [291 




the same interpretation and meaning, hold for more general, nonlinear Hamiltonian vector fields in the 
neighborhood of a saddle center-. . .-center equilibrium point, but obtaining the explicit formulae will require 
a normal form calculation for a given system. It is in this sense that restricting ourselves to quadratic 
Hamiltonians is, indeed, "without loss of generality", and it provides the method and approach for visualizing 
the same phase space structures governing reaction dynamics in more general systems. 



4 A model of the n-DoF system in the space of the integrals of 
motion 

For a system with more than 2 degrees of freedom, the dimensionality of the phase space structures discussed 
in Sec. [2] is too high to allow for a similar explicit visualization like for the 2-DoFmodels in Sec. [Sj Instead 
of projecting the phase space structures to the planes of the coordinates and conjugate momenta of the 



normal form like in Sec. 2.2 it is useful to present the phase space structures for n DoFin the space of the 



integrals /, J2, . . . , Jn- This can be accomplished using the momentum map M. defined in (12) in Sec. 2.3 



The discussion here follow s Waalkens et al.(2QQ8)| . 

As mentioned in Sec. 2.3 a fibre (the joint level set of the integrals /, J2,..., Jn in phase space) is 
called singular if it contains one (or more) irregular point (s). The image of all the singular fibres under the 
momentum map is called the bifurcation diagram. It is easy to see that the bifurcation diagram consists of 
the set of (/, J2, . . . , Jn), where one or more of the integrals vanish. In Fig. 34 we show the image of the 
energy surface with energy h > under the momentum map M. in the space of the integrals for n = 3 degrees 
of freedom. 




Figure 34: Sketch of the image of the energy surface of energy h > under the momentum map Ai in 
Equation ([l2| in the space of the integrals / and J/e, /c = 2, . . . , n, for the case of n = 3 degrees of freedom. 
The green/dark pink brindled piece of the image of the energy surface has / > 0; the turquoise piece has 
/ < 0. The intersections with the planes / = 0, J2 = and J3 = (pieces of which are visualised by 
semitransparent planes for clarity) form the bifurcation diagram of the energy surface. The image of the 
energy surface is not bounded in the direction of negative / as indicated by the dashed line at the bottom. 
The topology of the fibres A^~^(/, J2, J3) is indicated for the various points (/, J2, J3) marked by dots. The 
fibre of a point (/, J2, J3) with / 7^ consists of two disconnected manifolds as indicated by the factor of 2. 
The fibre of a point (/, J2, J3) with 7 = consists of a single connected manifold. 



The bifurcation diagram (of the energy surface) consists of the intersections of the image of the energy 
surface (the turquoise and green/dark ping brindled surface in Fig. 34) with one of the planes / = 0, J2 = 
or J3 = 0. Upon approaching one of the edges that have J2 = or J3 = the circle in the plane (^2,^2) or 
(^3,^3), respectively, shrinks to a point, and accordingly the regular fibres x R which form the Lagrange 
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manifolds A^j and Aj (see (13) and (14)) reduce to cylinders or 'tubes' S"*^ x R. At the top 

corner in Fig. 34 both J2 and J3 are zero. Here both circles in the centre planes (^2,^2) and (^'3,^3) have 
shrinked to points. The corresponding singular fibre consists of two lines, R, which are the forward and 
backward reaction paths, respectively (see also Fig. [4|. 

The fibres mentioned so far all have / 7^ and each consist of a pair of two disconnected components. 
For / < 0, one member of each pair is located on the reactants side and the other on the products side of 
the dividing surface. For / > 0, one member of each pair consists of trajectories evolving from reactants to 
products and the other member consists of trajectories that evolve from products to reactants. In fact the 
two members of a fibre which has / > are contained in the energy surface volume enclosed by the forward 
and backward reactive spherical cylinders Wf(h) and Wi)(h)^ see Fig. |4] For this reason we marked the piece 
of the image of the energy surface under the momentum map which has / > by the same green/dark pink 
colour in Fig. 34 that we used Fig. [4j Green corresponds to forward reactive trajectories and dark pink 
corresponds to backward reactive trajectories. Under the momentum map these trajectories have the same 
image. 

The light blue line in Fig. 34 which has / = is the image of the NHIM under the momentum map. For 
three degrees of freedom, the NHIM is a 3-dimensional sphere, and as we will discuss in more detail in Sec. [5] 
and indicated in Fig. 34 it is foliated by a one-parameter family of invariant 2-tori which shrink to periodic 
orbits, i.e. circles, S^, at the end points of the parameterisation interval. 



5 structure of the NHIM for n-DoF 

We have described how the NHIM, its stable and unstable manifolds, and the dividing surface, sit within the 
phase space and govern the reaction dynamics in several 2DoF geometrical models. However, we have not 
given any indication of the structure of the NHIM itself, except in the 2-DoF case, in which the NHIM is 
a single periodic orbit. In this section we discuss various ways of visualizing the geometrical structure and 
dynamics associated with the NHIM. 



Recall from Sec. 2.2 that in terms of the local decoupling provided by the normal form into reaction coor- 
dinates and bath modes, points on the NHIM, SNmM(^)' contain zero 'energy' in the reactive (saddle) mode 
(i.e., / = 0); all 'energy' is contained in the bath (centre) modes, which results in oscillatory, quasiperiodic 
motion with n — 1 independent frequencies. As also mentioned in Sec. |2.2[ this leads to the interpretation of 
the NHIM as the energy surface of an invariant subsystem with one degree of freedom less than the full system 
(a supermolecule located for a frozen reaction coordinate (and conjugate momentum) between reactants and 
products). 

Since, typically (for more than 1 DoF), there is more than one bath mode present, the total energy 
can be distributed in a variety of ways amongst the bath modes. For example, putting all of the energy 
into any one of the bath modes results in a single periodic orbit (often referred to as a Lyapunov orbit). 
More generally, the supermolecule of fixed energy h has energy distributed between all of the n — 1 bath 
modes. Since the integrals of motion are constant on trajectories, the amount of energy in each particular 
bath mode remains fixed during the motion. Such a molecule therefore undergoes quasiperiodic motion; it 
oscillates independently in the n — 1 different bath modes. The motion of a typical trajectory in the NHIM is 
quasiperiodic on an (n — l)-torus. Heuristically, one can think of the NHIM as the invariant surface made up 
of trajectories with all possible energy distributions between the bath modes (and no energy in the reactive 
mode). 

5.1 The NHIM (n = 3 DoF) 

As noted above, the NHIM for 2-DoF systems is simply a periodic orbit. This is an immediate consequence 
of the fact that there is only a single bath (or centre) mode. In this section we consider the next most 
complex case in some detail, the 3-DoF situation (2 bath modes), for which the NHIM is diffeomorphic to 
S^. However, first we will describe a way in which can be visualized in lower dimensions. This is helpful 
from the point of view of visual intuition. Our discussion in this section owes a great deal to the nice report 
of |ChishQlm(2QQQ)] . 
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Geometrical Construction of in Three Dimensions. To begin with, first note that in the usual 
manner, to draw a 1-sphere (i.e. a circle), we need two dimensons, to draw a 2-sphere we need three 
dimensions, and to draw a 3-sphere we need four dimensions. We will describe a method of drawing the 
3-sphere which requires only three dimensions. This particular representation will enable us to make a clear 
picture of the dynamics on the NHIM. The method can be viewed as a generalization of the presentation 
of a two-dimensional sphere, S^, in terms of two he misp heres (topological 2-balls, B'^) and an equator (a 
onedimensional sphere, S^) that we discussed in Sec. 3.1 For S^, one can form a similar representation by 



using two open 3-balls, 5^, to represent "Northern" and "Southern" hemispheres, and identifying points on 
their boundaries, S^, to make an "equator"; 



r3 

^North 



U ^Equator ^ ^South' 



(53) 

We win choose the "North" pole (N) to be at the centre of the North ball and the "South" (S) pole at the 



centre of the South ball. This model is illustrated in Fig. 35 



B 



North 



B 



South 




Figure 35: A model of the sphere formed by taking two open 3-balls, ^North ^South' identifying 
points on their boundaries; the "equator" ^1^^^^^^^. The arrow in illustrates how one passes from the 
chosen "North" to the chosen "South" pole at the center of each 3-ball. 



Partitioning 'Energy' Between the Two Bath Modes Recall that (for 3 DoF) there are 2 bath modes. 
For a given trajectory in the NHIM, all of the 'energy' could be in the first bath mode (resulting in a periodic 
orbit) or in the second bath mode (resulting in another periodic orbit), or the energy could be distributed, 
or partitioned, between the two bath modes. 

'a) we show the two open balls that comprise our model of as described above. We choose 



In Fig. 
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a trajectory having the property that all of the energy is in the first bath mode, i.e., J3 = and J2 = JJ^^^, 
where J^^^^ is defined by the energy equation, i^^NF(0, J2^^^,0) = h > (see Sec. 2.2). This corresponds to a 
periodic orbit that moves along the "vertical" circle that passes through the centers of the two balls. As the 
periodic orbit exits the North ball at the top, it enters the South ball from the top, exiting the South ball at 
the bottom, and entering the North ball at the bottom. 

Suppose now that we change this distribution of energy very slightly, by taking a small amount of energy 
out of the first bath mode and putting it into the second, which we illustrate in Fig. [36]^b) (and where we 
leave it to the reader to consider how a trajectory leaves one ball through its surface and enters the other 
ball). The motion still stays close to the original circle, but there is now a second circular (oscillatory) motion 
"around the periodic orbit". Thus the motion now takes place on a thin 2-torus. Any initial condition on 
the 2-torus corresponds to a quasiperiodic trajectory that remains on the 2-torus. 

We continue this procedure of taking energy from the first bath mode and putting it into the second bath 
mode, until all of the energy is in the second bath mode (Fig. 36 ^e)), which results in a periodic orbit moving 
in a "horizontal" circle along the surface of the two balls. 
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(e) 



Figure 36: Illustration of the invariant tori in for n = 3 degrees of freedom, (a) representation 

of the periodic orbit corresponding to local integrals of motion (/, J2, J3) = (0, J2^^^,0), shown in both 
hemispheres of the model, (b) a nearby torus with J3 = e with < e <C JJ^^^, (c) the case J2 = J2^^^/2, (d) 
a torus with J3 = J^^^ — e which is close to the periodic orbit corresponding to (/, J2, J3) = (0,0, Jf^^^), 
shown in (e). 
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In this way we see that is foliated by a one-parameter family of 2-tori. Trajectories on the 2-tori 
correspond to quasiperiodic motion with energy contained in both bath modes. The "extreme" cases corre- 
sponding to all the energy contained in one of the bath modes corresponds to periodic orbits (the "Lyapunov 
orbits"). We illustrate this family of 2-tori in Fig. 37 




Identify 



Figure 37: The family of invariant 2-tori that fill the NHIM, S^j^jj^, illustrated according to our method 
for visualizing in three dimensions. One can think of each torus as being cut in half, with one half in 
the North ball and the other half in the South ball, and the "edges" of the two halves touching touching 
the surface of the 3-ball, where they are identified. In this way, a quasiperiodic trajectory winding around a 
2-torus continually moves from one 3-ball to another as it encounters the edge of the torus where the torus 
is cut into two pieces, which are on the surfaces of the 3-balls (i.e. the "equator" of S^him)- 



5.2 The NHIM, SnhimI^)' forn > 3 DoF 

The method for visualizing the NHIM for 3 DoF described in the previous section is not as useful for n > 3 
DoF. Nevertheless, the dynamics on the NHIM have the same characteristics-quasiperiodic trajectories with 
the energy partitioned amongst the n — 1 bath modes. Hence, as mentioned in Sec. [4j for the general DoF 
case, a description of the NHIM can be given in terms of the geoemtry of the space of centre actions (or 
"bath mode" actions) J2, • • • , -^n- This is what we will now describe. 

Classically, the space formed by the centre actions is the (n — l)-dimensional positive cone; equivalently, 
the positive orthant (the generalization of a quadrant in the plane or an octant in three dimensions to 
arbitrary dimensions) given by J2, Jn > in R^~^. In general, the NHIM is represented in this space by 
a (topological) (n — 2)-simplex having (n — 1) vertices, one for each centre action. The vertices are always 
located on the coordinate axes; each vertex thus corresponds to (n — 2) of the integrals being equal to 
zero. Note that the projection of the NHIM to this space of actions will generally be a nonlinear embedding 
of a simplex; the embedding is only linear (that is, all cells of the simplex are flat) in the case where the 
vector field defined by Hamilton's equations is linear. We will denote this embedded simplex by 

A=[2,3,...,n], (54) 

using the index of each centre action to label a vertex. Within this simplex, each {k — l)-cell (for every 
k = 1, . . . , (n — 1)), denoted [ji, . . . , j/e], represents a (/c — l)-parameter family of invariant /c-tori which we 
denote by 

(55) 

Here ji, , . . . , j/c are mutually different elements of {2, . . . , n} which indicate the center actions J^^ , . . . , J^^ 
are different from zero. The family is parameterised, for example, by the centre actions J^^, . . . , with 
the value of Jj^ being determined by the energy equation, as we will describe below. 



38 



For example, each vertex (i.e., O-ceh), denoted [j] (for j = 2, ...,n), represents a single periodic orbit 
Tj = S^, in which all actions except Jj are set to zero and, therefore, the remaining action is fixed at a value 
Jj — Jj^^^, which solves the following equation that corresponds to restriction to a single energy surface, 

Knf(/,J2,..., Jn) =/^>0, (56) 

subject to setting = for all ^ = 2, . . . , n, ^ 7^ j. The periodic orbit, Tj, projects to a single circle of 
radius y^J^^ in the ((7/e,p/c)-plane of the normal form coordinates and to the origin in all other normal form 
planes, {qi.Pi) for i^^j. 

Each edge [^1,^2] (for ji 7^ 32) of the simplex represents a 1-parameter family of invariant 2-tori, in which 
all actions except Jj^ and Jj^ have been set to zero. Without loss of generality, we can choose Jj^ as a 
parameter; each value Jj^ G [0, Jj^^^] fixes the value of Jj^ by virtue of the above iso-energetic condition 
subject to the modified restriction that 

J^ = 0for€^ji,i2. (57) 



We denote this 1-parameter family of invariant 2-tori in Sjs^^jj^(/i) by T^^ 



By allowing k centre actions to be non-zero in the NHIM, corresponding to a (/c — l)-cell [ji, ^2 • • • , j/c] of 
the simplex, we can, without loss of generality, take Jj^, Jj2, . . . , Jjk_i as parameters which fix the value of 



Jj^ via the condition 



i7NF(/,J2,...,Jn) = /^>0, (58) 

subject to setting 



This defines a (A: — 1) -parameter family of invariant /c-tori, denoted T^^ 



= for ^ 7^ ji, j2, . . . , jfe. (59) 

jlj2---jk ' 



A schematic representation of the projection of S^him(^) ^^^^ space of centre actions J2, -^3, • • • , -^n 



for n-DoF is shown in Fig. 38 for (i) n = 3, (ii) n = 4, (iii) n = 5. In (iii) we show a schematic 3-dimensional 



projection of the 4-dimensional space of centre actions. 



Directional flux across the dividing surface. In [Waal kens fc Wiggins (2004)] it is shown that the 
(directional) fiux across a given iso-energetic dividing surface {h) (that is, the equal and opposite fiuxes 

through the forward and backward dividing surface hemispheres B'^~^'^{h) and B'^^{h)) may be reduced to 
a certain integral taken over S^^^(/i). In the representation of the NHIM in in terms of the centre integrals 
given above, this fiux corresponds to the (n — l)-dimensional volume of the finite region of the positive 
orthant bounded by the (n — 2)-simplex representing the NHIM times {2t:)^~^ . That is up to a factor, the 
fiux corresponds to the volume of an (n — l)-simplex formed by adding the origin to the vertex set of the 
NHIM in the space of integrals. In Fig. |38]for (i) and (ii), the (n — l)-dimensional volume F corresponding 
to the magnitude of the directional fiux across the dividing surface hemispheres B' ^~^'^{h) and B'^^(h) i s 
indicated. Applications of this in both the classical and quantum case can be found in [Waalkens et al.(2QQ8)| . 

6 Conclusions 

This paper is concerned with strategies and approaches for visualizing the phase space structures that govern 
reaction dynamics. In Section |2] we reviewed recent results that develop the phase space geometry near 
an equilibrium point of Hamilton's equations of saddle x centre x • • • x centre stability type ('saddle' for 
short) in n degree-of-freedom Hamiltonian systems. The key technique used was the method of Poincare- 
Birkhoff normal forms. This method provided a transformation (and its inverse), valid in a phase space 
neighborhood of the equilibrium point, to normal form coordinates in which the phase space structures could 
be explicitly realized (and then mapped back to the original coordinates via the inverse of the normal form 
transformation). The key phase space structure was seen to be a normally hyperbolic invariant manifold 
(NHIM), which is the higher dimensional analogue of a "saddle point". The NHIM is the anchor upon 
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h h 

(i) n = 3 degrees of freedom. 

J3 




4 J4 

(ii) n = 4 degrees of freedom. 

J3 

[3]ft. [3,4,5] 
[3.4] 




(iii) n = 5 degrees of freedom. 

Figure 38: Schematic representation of the projection of S^^m(^) ^^^^ space of centre actions, 
J2, J3, . . . , Jn, for n DoF giving a nonhnear embedding of a simplex with (n — 1) vertices, [j], j = 2, . . . , n, 
which encloses a volume proportional to the directional flux, F. 



which a dividing surface having the "no-recrossing" property and minimal flux is constructed. Moreover, the 
stable and unstable manifolds of the NHIM are the "conduits" for forward and backward reactive trajectories 
(trajectories evolving from reactants to products, and vice- versa). More precisely, one the two branches of 
the stable manifold can be glued together with one of the two branches of the unstable manifold to form a 
cylinder (denoted the forward reactive cylinder in this paper) that enclose all forward reactive trajectories; 
and the other branches of the stable and unstable manifolds can be glued together to form a cylinder (the 
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backward reactive cylinder) which encloses ah backward reactive trajectories. In the case where a generic 
non-resonance condition holds on the pure imaginary eigenvalues associated with the center directions it 
follows that in the normal form coordinates the Hamiltonian (when the series representing the normal form 
Hamiltonian is truncated at some finite order) has n integrals of the motion and thus is integrable. This 
integrable structure allows for a natural decoupling of the dynamics into a single reactive mode and n — 1 bath 
modes. Moreover, the integrable structure provides a foliation of the NHIM by a family of lower dimensional 
tori and a foliation of a neighborhood of the equilibrium point into Lagrangian cylinders. 

The "decoupling" of the vector field and the associated dynamics masks the complete picture of reaction 
governed by the phase space structures in the energy surface. In Section |3] we provide three different 2 
degree-of-freedom models that allow us to visualize the energy surface, the NHIM and its stable and unstable 
manifolds, the dividing surface, and the Lagrangian cylinders. These give a complete picture of how the 
phase space structures govern the dynamics of trajectories near a saddle. By contrast, the role of the phase 
space structures is typically obscured in their common representation as projections to configuration. For 
example, the projection of the volume enclosed by the forward and backward cylinders might not lie inside of 
the projection of the forward and backward cylinders, and as a consequence forward and backward reactive 
trajectories might seemingly leave the forward and backward reactive cylinders (for an example, see the study 
of the planar Hill problem in [Waalkens et al.(2QQ5b)j ). 

For the general case of n > 2 degrees of freedom systems, we use the integrals alone to provide a model 
for visualizing the phase space structures (see Section [4] ) . This model does not contain explicit information 
on trajectories, but it can be easily related to the behavior of trajectories by utilizing the foliation of the 
neighborhood of the NHIM by Lagrangian cylinders. 

We moreover discussed in some detail the foliation of the NHIM for 3 degrees of freedom. Here the NHIM 
is a threedimensional sphere, that is foliated by a one-parameter family of invariant 2-tori. We visualized 
this foliation using the 3-ball model of (see Section [5|. For the general case of n > 3 degrees of freedom, 
we presented the NHIM as a simplex in the space of the integrals. This is of particular importance, since, 
up to a prefactor, the volume enclosed by the simplex in the positive orthant of the integral space gives the 
directional fiux. 

We conclude by emphasizing that we focused the discussion of the phase space structures and their 
visualization to the neighborhood of the saddle. In fact, depending on the global topology of the energy 
surface, the local partitioning of the energy surface by the dividing surface constructed from the NHIM into 
a reactants and a products component might not extend to a global partitioning. For example, there could 
be more than one saddle point which controls the access to and escape from a certain phase space region 
which requires a more careful definition of a reactants and a products region. Moreover, in this paper we 
did not discuss the global geometry of the stable and unstable manifolds of the NHIM. In fact the stable 
and unstable manifolds might extend to regions far away from the neighborhood of the saddle. Even if there 
is only a single saddle point, the global structure might be very complicated since the stable and unstable 
manifolds might wrap around in a complicated manner and intersect each other to form a highdimensional 
homoclinic tangle. If there are several saddles the stable and unstable manifolds of NHIMs of different saddles 
might intersect. Since the stable and unstable manifolds are of codimension 1 in the energy surface they still 
act as impenetrable barriers globally, and also control the global dynamics in a possibly complicated way. 
The details of the dynamics then strongly depends on the system under consideration. For a few examples, 
we refer to [Waalkens et al.(2QQ4a)[ [Waalkens et al.(2QQ4b)[ [Waalkens et al.(2QQ5a)l [Waalkens et al.(2QQ5b)| 



Waalkens et al.(2QQ5c)| 
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